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UNSTEADY FLOW METER 

This present invention concerns the area of measurement 
of the instantaneous flow of a fluid in unsteady motion. 
5 It finds particular but not limiting advantageous 

application in the areas of process engineering and the motor- 
vehicle industry. 

More particularly, it finds advantageous application in 
the control and regulation of internal-combustion engines, 
10 test centres and research laboratories, and all the current 
applications of flow meters from previous art. 

More particularly, it finds advantageous application in 
all the situations where the flow is unsteady, and knowledge 
of the flow requires the integration of speed in space and 
15 time. 

GENERAL PRESENTATION OF PREVIOUS ART 

Flow measurement is an essential element in the 
optimisation of industrial processes. 

20 At the present time, in order to ascertain the flow of a 

fluid in unsteady motion at a given instant, it is generally 
necessary to employ measurement techniques that enable the 
speed of this fluid to be determined. Once a set of speeds 
have been determined for different instants in time, this set 

25 of speeds is then integrated spatially in order to ascertain 
changes to the flow over time. 

Many systems and methods have already been proposed as 
the means to obtain the flow of a fluid in unsteady motion by 
determining the speed of this fluid. 

30 By way of an example, one can mention the LDV and PIV 

methods, and the systems which are associated with them. 

The PIV method, or Particle Image Velocimetry, is based 
upon measuring the movement of small particles. Determining 
the speed of the fluid at a given instant is effected in the 

35 following manner. The moving fluid is seeded with particles, 
and is then illuminated by a laser. The laser emits two pulses 
offset by a short time interval. A CCD camera synchronised to 
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the laser records an image for each pulse from the laser. The 
two images obtained are then processed. The movement of the 
particles between the first and the second image is measured 
using a correlation technique. By dividing the movement of the 
5 particles between the two images by the interval separating 
the recording of each of these images, one obtains the speed 
of the moving fluid. This determination of the speed of the 
fluid is then used to deduce the flow of this fluid. 

The LDV method, or Laser Doppler Velocimetry, is based 

10 upon the measurement of a Doppler frequency. The moving fluid 
is seeded with particles and is then illuminated by means of 
two lasers whose beams cross. The dispersed light coming from 
the particles passing through the intersection of the two 
laser beams then produce a Doppler frequency which is 

15 proportional to the speed of the particle. Calculation of this 
speed enables the flow of the moving fluid to be found. 
Viewing of the flow is achieved in the same manner as in the 
PIV case. 

These methods and systems have the following drawbacks 
20 however: 

the PIV and LDV methods cannot be used to determine 
the flow of a fluid in real time. In fact, these techniques 
require post-processing of the acquired images in order to 
obtain the speed of the fluid. 
25 - The use of such techniques is onerous and expensive. 

These methods are also poorly resolved in time, and 
therefore have a small passband (generally less than 10 Hz for 
the PIV method) . 

the fluid has to be seeded with particles. 
30 - the fluid and the conduit must be transparent. 

Another method used to obtain the flow of a moving 
unsteady fluid from its speed employs the Pitot tube. 

The Pitot tube is an instrument that is used to determine 
the speed of a moving fluid by measuring a pressure 
35 difference. In order to determine the speed of a moving fluid 
in a conduit, it is necessary to ascertain the direction of 
the flow. The Pitot tube is inserted into the conduit so that 

WASH_1667796.1 



the tube is parallel to the direction of the flow, and so that 
the front of the Pitot tube is facing into the flow. A first 
pressure is obtained by taking a first pressure reading at the 
lateral wall of the Pitot tube. This first measured pressure 
5 is a static pressure. A second pressure is obtained by reading 
the pressure at the front of the tube. This pressure is the 
total pressure of the fluid. Measuring the difference between 
these two pressures enables the speed of the moving fluid to 
be calculated. From the calculation of this speed, it is 
10 possible to obtain the flow of the moving fluid. 

This method and this system have the following drawbacks 
however : 

the direction of the flow must be known and constant 

in order to position the Pitot tube. 
15 - when the direction of the flow varies, the total 

pressure measurement is adversely affected, and this causes an 

error in the acquisition of the speed measurement. 

- this type of system is intrusive, and therefore 

disrupts the flow. 
20 A final method, used to ascertain the flow of a fluid in 

unsteady motion from its speed, is the so-called hot-wire 

method . 

The hot-wire method is based upon heat transfer around a 
wire heated by an electrical circuit. One places into the 

25 moving fluid a wire that is raised by eddy-current effect to a 
temperature which is greater than the temperature of this 
fluid. There then occurs an exchange of heat by convection, 
where the wire is cooled by the moving fluid. It is then 
possible to determine the speed of the moving fluid by 

30 calculating the power that has to be supplied to the wire in 
order to keep it at a constant temperature. When the speed of 
the fluid is known, it is then possible to deduce the flow of 
the fluid in unsteady motion. 

This method and this system have the following drawbacks 

35 however: 

the method is not automated and is insensitive to 
the direction of fluid flow. 
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measurements in the liquids, which are conducting 
most of the time, require that the heated wire (heated by the 
electrical circuit) should be isolated electrically, 
the system has to be calibrated. 
5 - the wire ages quite rapidly and this calls for 

costly maintenance. 

the deployment, and in particular the dimensions, of 
the wire results in serious fragility of the system. 

it is an intrusive method, and therefore liable to 
10 disrupt the flow. 

One objective of this present invention is to provide a 
method and a system for measuring, in real time, the movement 
of a fluid in unsteady motion, and which overcomes most of the 
aforementioned drawbacks . 

15 

PRESENTATION OF THE INVENTION 

The invention concerns a measurement system for the 
measurement, in real time, of the instantaneous flow of a 
fluid in steady or unsteady motion in a conduit, characterised 
20 in that it includes: 

a flow-velocity measuring device placed in the 
conduit, where the said flow-velocity measuring device is 
equipped with two pressure take-offs in its wall; 

a means of measuring a pressure difference, designed 
25 to be connected to the two pressure take-offs; 

a calculation resource designed to calculate the 
flow in real time by solving an equation that relates the 
instantaneous flow to the pressure difference, with the latter 
being positive or negative in the said equation depending on 
30 variations in the speed of fluid flow in the conduit and/or 
the direction of the fluid flow. 

As will be explained later, the flow measurement system 
of this present invention allows the industrial monitoring of 
even highly pulsed flows, combing low price, simplicity of 
35 implementation and reliability. This flow measurement is 
performed in real time and takes account of the direction of 
the flow. In addition, the flow measurement system of this 
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present invention is non-intrusive, and does not require 
calibration . 

Preferred, but not limiting, aspects of the flow 
measurement system according to the invention are as follows: 
5 - the equation includes a term representing the 

pressure difference, a term representing the instantaneous 
flow, and a term representing the differential coefficient 
with time of the instantaneous flow, where each of these three 
terms can be positive or negative. 
10 - the equation taking account of the direction of the 

flow is a differential equation of the form: 

dq(t)/dt + a(q(t)) = P x Ap(t), 

15 - where q(t) is the wanted instantaneous flow, 

dq(t)/dt is the differential coefficient with time of the 
wanted instantaneous flow, <x(q(t)) is a function that depends 
on the geometry of the system, the fluid and the flow q(t), (3 
is a coefficient that depends on the geometry of the device, 

20 and Ap(t) is the measured instantaneous pressure difference. 

the system also includes a temperature measuring 

probe . 

the calculation resource is designed to ascertain 
the density of the fluid by having the temperature measured by 
25 the temperature measuring probe, and calculating the 
instantaneous mass flow of the fluid. 

the system also includes a probe for measuring the 
absolute static pressure. 

the calculation resource is designed to calculate, 
30 in real time, the instantaneous mass flow of a fluid that is 
compressible in real time, by measuring the absolute static 
pressure, by measuring the temperature, and by solving the 
equation relating the instantaneous flow to the pressure 
difference, where the said equation takes account of the 
35 direction of the fluid flow. 

the flow-velocity measuring device is a narrowed 

tube . 
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the flow-velocity measuring device is a diaphragm, 
the flow-velocity measuring device is a venturi. 
the means of measuring the pressure difference is a 
differential pressure sensor connected to the two pressure 
5 take-offs. 

the means of measuring the pressure difference is a 
set of two relative pressure sensors connected to the two 
pressure take-offs . 

the calculation resource is an electronic 
10 calculator. 

the electronic calculator used to calculate the 
instantaneous flow in real time, irrespective of the 
fluctuations in the flow, is an analogue or digital electronic 
calculator . 

15 - the calculation resource includes a first amplifier 

connected to a first input of a subtractor, an integrator 
connected to an output of the subtractor, a feedback loop 
connected between an output of the integrator and a second 
input of the subtractor, where the feedback loop includes a 

20 module to produce the absolute value function, connected to 
the output of the integrator, and a multiplier connected by a 
first input to an output of the module and by a second input 
to the output of the integrator. 

This present invention also concerns a process for 

25 measuring the flow of a fluid in steady or unsteady motion in 
a conduit, characterised in that it includes the following 
stages: 

the measurement, irrespective of fluctuations in the 
flow, of a difference between two pressures. 

30 - calculation of the flow of the fluid by solving an 

equation that relates the flow and the measured' pressure 
difference, where the latter is either positive or negative in 
the said equation depending on variations in the speed of 
fluid flow in the conduit and/or the direction of the fluid 

35 flow. 

Preferred but non limiting aspects of the flow 
measurement process according to the invention are as follows: 
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the second stage of the process consists of 
calculating the flow of the fluid in unsteady motion by 
solving a differential equation relating the flow and the 
measured pressure difference : 

5 

dq(t)/dt + a(q(t)) = p x Ap(t), 

where q(t) is the wanted instantaneous flow, 
dq(t)/dt is the differential coefficient with time of the 
10 wanted instantaneous flow, a(q(t)) is a function that depends 
on the geometry of the system, the fluid and the instantaneous 
flow q(t), P is a coefficient that depends on the geometry of 
the device, and Ap(t) is the measured instantaneous pressure 
difference . 

15 - the direction of the flow is included in the term 

<x(q(t)) which is dependent on the geometry of the system and 
the flow q ( t ) . 

the process also includes a stage for acquisition of 
the fluid temperature. 

20 - the process also includes a stage for acquisition of 

the absolute static pressure of the fluid. 

the temperature of the fluid and the absolute static 
pressure of the fluid are included in the equation taking 
account of the direction of the flow. 

25 - the process allows the calculation, in real time, of 

the instantaneous flow of a compressible fluid in unsteady 
motion in a conduit, where the said calculation of flow is 
achieved by solving an equation that relates the flow to the 
pressure difference, the absolute static pressure of the fluid 

30 and the fluid temperature. 

PRESENTATION OF THE FIGURES 

Other characteristics, objectives and advantages of this 
present invention will emerge from the description that 
35 follows, which is purely illustrative and non limiting, and 
which must be read with reference to the appended drawings, in 
which : 
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Figure 1 illustrates a view in section of one method 

of implementation of this present invention; 

Figure 2 illustrates a block diagram to be 

implemented in order to determine the flow of a fluid 
5 according to the process of this present invention; 

Figure 3 illustrates an example of a numerical 

algorithm to be executed in order to determine the flow of a 

fluid according to the process of this present inventions- 
Figure 4 illustrates a comparison between flow 
10 measurement by the process of this present invention and flow 

measurement by the hot-wire method. 

DESCRIPTION OF THE INVENTION 

One objective of this present invention is to establish a 

15 process and a system allowing measurement of the flow of a 
fluid in unsteady motion in a conduit, even when the direction 
of the flow varies over time, and in which flow measurement 
takes account of the direction of the flow. 

To this end, the inventors have established a relation 

20 between the instantaneous flow and the pressure of a fluid in 
unsteady motion in a conduit, which enables a flow to be 
calculated in real time, and where the value of this flow can 
be positive or negative, so that the sign of the calculated 
flow indicates the direction of the fluid flow. 

25 With this in mind, the inventors started with the 

aerodynamic equations. By firstly applying the usual 
approximations for flows in a conduit to the equation for 
total kinetic energy, and then replacing each speed-of-f low 
term by its expression as a function of the flow in the total 

30 kinetic energy equation, the inventors arrived at an equation 
relating the instantaneous flow and the pressure for a fluid 
in semi-unidirectional unsteady motion in a conduit. 

This equation is a differential equation of the following 

form: 

35 dq(t)/dt + <x(q(t)) = (3 * Ap(t), 

where : 

q(t) represents the wanted flow by volume, 
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dq(t)/dt represents the differential coefficient with 

time of the wanted flow, 
<x(q(t)) represents a function that depends on the 
geometry of the system, the fluid and the flow 
5 q(t), 

P represents a coefficient related to the 

geometry of the device, 
Ap(t) represents the measured pressure difference. 

The fact that the equation described above is a 
10 differential equation enables one to ascertain the flow by 
volume in real time, and this could not be done correctly with 
an algebraic equation. 

This differential equation comes directly from the fluid 
mechanics equations. This equation has the particular feature 
15 that it takes account of the direction of the flow. The second 
term on the left side, a(q(t)), is a function of the flow 
configured by the geometry of the device and by the direction 
of the flow. This property allows the flow to be measured even 
when the flow direction reverses. 
20 The presence in the equation of terms whose sign varies 

according to variations in the speed of fluid flow in the 
conduit and/or the direction of the flow (term Ap(t)), and of 
terms whose sign varies according to the direction of the 
fluid flow in the conduit (terms dq(t)/dt and a(q(t))) enables 
25 one to calculate a positive or negative flow, whose sign 
indicates the direction of the fluid flow. 

In particular, in the equation described above, the term 
Ap(t) is employed just as it is, with no absolute value. The 
sign attributed to the pressure difference, Ap(t), varies in 
30 particular as a function of the variation of speed of fluid 
flow and/or according to the direction of the fluid flow. The 
pressure difference Ap(t) is therefore positive or negative, 
as for example when the fluid flowing in the conduit 
accelerates or decelerates, and/or when the fluid flows in the 
35 conduit in one direction or in the opposite direction. 
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In addition, in the equation described above, the terms 
dq(t)/dt and ot(q(t)) are also used as they are, with no 
absolute value. 

This enables one to calculate a flow whose value can be 
5 positive or negative, in contrast to the devices of previous 
art, which are capable of calculating only positive flows and 
not negative flows. 

The sign of the calculated flow enables account to be 
taken of the direction of the flow, and the inversions of sign 
10 of the calculated flow allow account to be taken of the 
reversals of direction of the fluid flow. . 

For example, if the calculated flow, at t = ti, is equal 
to (-20 m 3 /s) , this means that the fluid flows from left to 
right in the conduit. If the calculated flow at t = t 2 is 
15 equal to (+5 m 3 /s) , then this means that there has been a 
reversal of the direction of the flow between t = ti and 
t = t 2 , and that the fluid is flowing from right to left in 
the conduit at t = t 2 . 

The reader will have understood that "absolute value" 
20 refers to the square root of the square of a number, where 
this value can be positive or zero but never negative. 

This differential equation has an original form which has 
the property of always converging on condition that the 
initial condition is correctly chosen in the order of 
25 magnitude of the flow to be found. 

Preferably, the initial condition will be chosen at the 
moment when the flow changes direction, so that q(t=0) = 0. 

DESCRIPTION OF ONE METHOD OF IMPLEMENTATION OF THE INVENTION 

30 As shown in figure 1, the inventors have created a system 

which is designed to measure, in real time, the flow of a 
fluid in a conduit on the basis of a pressure difference. This 
system is a flow meter, and is composed mainly of a flow- 
velocity measuring device 1 inserted into the conduit, a 

35 pressure-difference measuring resource 2, and a calculation 
resource 3 designed to calculate the flow in real time. 
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The flow-velocity measuring device 1 illustrated in 
figure 1 is a narrowing tube. It is a profiled conduit with a 
cylindrical first part 4 that has the section SI. This first 
part 4 ends in a second convergent part 5 which leads to a 
5 third cylindrical part 6 of smaller section than SI. Two 
static pressure take-offs A and B are placed respectively on 
the wall of the cylindrical first part 4 and on the wall of 
the third cylindrical part 6 of smaller section. 

The pressure-difference measuring resource 2 illustrated 

10 in figure 1 is a differential pressure sensor. The pressure- 
difference measuring resource is connected to the two static 
pressure take-offs A and B of the flow-velocity measuring 
device 1. The pressure-difference measuring resource 2 enables 
the difference between the two static pressures to be 

15 measured. 

The calculation resource 3 is an electronic calculator. 
This calculation resource 3 is programmed to execute the 
algorithm for solving the equation to calculate the flow, and 
the sign of this flow takes account of the direction of the 

20 flow. The calculation resource is designed to calculate the 
flow in real time by solving an equation relating the 
instantaneous flow to the pressure difference, this pressure 
difference being positive or negative in the said equation, 
depending on variations in the speed of fluid flow in the 

25 conduit and/or the direction of the fluid flow. The equation 
relating the instantaneous flow to the pressure difference has 
the advantage of allowing the calculation resource to 
calculate a flow regardless of the direction of the current, 
to detect the direction of fluid flow, and to detect changes 

30 in the direction of the fluid flow. This calculation resource 
can be digital or analogue, as illustrated in figures 2 and 3. 

The principle of operation of the measuring system of 
figure 1 is as follows : 

The fluid flows in a conduit and passes through the flow- 

35 velocity measuring device 1. 

The pressure-difference measuring resource 2 measures the 
difference between the static pressure PI acquired at pressure 
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take-off A, and the static pressure P2 acquired at pressure 
take-off B. The pressure difference measured by the measuring 
resource 2 is transmitted in the form of an electrical signal 
to the input of the calculation resource 3. 
5 The calculation resource 3 calculates the flow of the 

fluid in real time from the pressure difference received as 
input. In fact, the calculation resource 3 is designed to 
solve the equation relating the flow to the pressure 
difference, where this equation allows calculation of the flow 

10 in real time and also gives the direction of the fluid flow. 

At the output of the calculation resource 3, one obtains 
the flow of the fluid in unsteady motion. 

The device of figure 1 therefore allows the measurement 
in real time of the flow of a fluid in a conduit. Measurement 

15 of this flow requires only two measurements of static 
pressure, and this is achieved by a single differential 
pressure sensor. Creating the device of this present invention 
is therefore simplified, and the presence of a single 
differential pressure sensor in the device results in a 

20 reduced cost of the said device. 

In addition, in contrast to the systems of the previous 
art, the device of this present invention allows reversals of 
direction of the flow, and takes account of the direction of 
the flow. 

25 In addition, the device of figure 1 does not require flow 

calibration, and the presence of only two pressure take-offs 
in the wall does not disrupt the flow, in contrast to a 
measurement system by hot wire or Pitot tube for example, 
where the probe is inserted directly into the conduit. 

30 The block diagram of figure 2 illustrates in principle 

the integration of the differential equation used to calculate 
the flow, and taking account of the direction of the flow, by 
means of a closed-loop circuit. This diagram comprises an 
integrator closed by a negative feedback loop which provides 

35 asymptotic stability of the integral form of the differential 
equation, that is of the equation to calculate the flow in 
real time and giving the direction of the flow. 
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The calculation resource 3 can be analogue. As 
illustrated in figure 2, the calculation resource is composed 
of a first amplifier 7, which can be purely passive, a 
subtractor 8, an integrator 9, and a feedback loop that 
5 includes a circuit 10 (or module) producing the absolute value 
function, a multiplier 11 with two inputs, and a second 
amplifier 12. The feedback loop enables transmission to the 
input of the subtractor 8 of a value that is a function of the 
calculated flow obtained at the output of the calculation 

10 resource (that is at the output of the integrator 9. 

The input to the first amplifier 7 is the input to the 
calculation resource. The first amplifier 7 amplifies the 
input signal of the calculation resource by coefficient p. 
This coefficient is the value of the coefficient P of the 

15 differential equation. This coefficient is a function of the 
geometry of the device, and more particularly of the geometry 
of the flow-velocity measuring device. The value of this 
coefficient is determined once only. This process can be 
effected in different ways, such as: 

20 - by comparative calibration for example, 

by directly evaluating the flow of a real or 
simulated fluid, constant or periodic, of period T (in the 
simulated case, this requires knowledge of the pressure 
calibration of the differential sensor) , 

25 - by calculation. 

The first amplifier 7 receives as input the input signal 
from the calculation resource. This input signal from the 
calculation resource is a signal representing the pressure 
difference. It is taken from the pressure-difference measuring 

30 resource 2. The output of this first amplifier 7 goes to the 
subtractor 8 . 

The subtractor 8 is used to obtain the difference between 
the output signal from the first amplifier 7 and the signal 
obtained from the feedback loop. This signal obtained from the 
35 feedback loop corresponds to the term <x(q(t)) of the 
differential equation (here we have <x(q(t)) = a * q(t) x 
lq(t) | ) . 
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The output from the subtractor 8 is used as input to the 
integrator 9. This integrator 9 is then use to integrate the 
output signal from the subtractor 8 in real time. The output 
signal from the subtractor 8 corresponds to the differential 
5 coefficient with time of the wanted flow (In fact from the 
equation we have dq(t)/dt = (3 x Ap(t) - <x(q(t))). The signal 
appearing at the output of the integrator 9 corresponds to the 
wanted flow (q) . This is the output from the calculation 
resource . 

10 The signal appearing at the output of the integrator 9 is 

sent to the input of the multiplier 11. It is also sent to the 
input of the circuit 10 producing the absolute value function. 
The output from the circuit 10 producing the absolute value 
function is used as the second input to the multiplier 11. 

15 The multiplier 11 effects the multiplication of the 

output signal from the integrator 9 by the output signal from 
the circuit 10 producing the absolute value function. The 
output of the multiplier ■ 11 is then used as input to the 
second amplifier 12. 

20 The second amplifier 12, which can be purely passive, 

amplifies the output signal by a coefficient a. This 
coefficient is a function of the geometry of the device and 
more particularly of the geometry of the fluid velocity 
measuring device 1. It is also a function of the 

25 characteristics of the fluid. The value of this coefficient is 
determined once only. This determination can be effected in 
different ways, such as, for' example: 
- by comparative calibration, 

by directly evaluating the flow in a real or 

30 simulated periodic flow of period T (in the simulated case, 
this requires knowledge of the pressure calibration of the 
differential sensor) , 

by calculation. 
In practice, determination of parameters a and P in 

35 operation requires adjustments which are interdependent. 

The principle of operation of the calculation resource of 
figure 2 is as follows: 
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As indicated previously, the differential equation 
relating the flow of a fluid in unsteady motion to a pressure 
difference is of an original form which has the property of 
always converging on condition that the initial condition is 
5 correctly chosen, of the order of magnitude of the mean flow 
(e.g. q(t=0) = 0 when the flow changes direction). This means 
that given an initial condition on the output signal from the 
integrator, this output signal will converge to a value 
representing the wanted flow. 

10 To begin with, there will be a zero signal, for example, 

at the output of the integrator (q(t=0) = 0). The calculation 
resource receives as input a signal representing the pressure 
difference between the two pressure take-offs. This signal is 
amplified by the first amplifier 7. The output signal from 

15 this first amplifier corresponds to the term p * Ap(t) of the 
differential equation . 

The output signal from the first amplifier 7 is 
transmitted to the first input of the subtractor 8. The signal 
obtained from the feedback loop arrives at the second input of 

20 the subtractor 8 (equal to zero to begin with, since the 
initial condition chosen is q(t=0) = 0). The subtractor 8 
effects the subtraction of these two signals. The output 
signal from the subtractor correspond to the term dq/dt in the 
differential equation . 

25 The output signal from the subtractor 8 is sent to the 

integrator 9. This signal is integrated . by the integrator 9 
which delivers as output a signal representing the wanted flow 
(q) . 

The output signal from the integrator 9 is sent to the 
30 output of the calculation resource and as input to the 
feedback loop. The input signal to the feedback loop is sent 
to one of the inputs of the multiplier 11, and also to the 
input of the circuit 10 producing the absolute value function. 
The output signal from the integrator 9 and the output signal 
35 from the circuit 10 producing the absolute value are 
multiplied by the multiplier 11. 
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The output of the multiplier 11 is then sent to the 
second amplifier 12 which amplifies it. The output of the 
second amplifier 12 corresponds to the term a(q(t)) of the 
differential equation . 
5 The output signal from the second amplifier is sent to 

the second input of the subtractor. This signal takes account 
of the direction of the flow . 

The subtractor effects the subtraction of the output 
signal from the first amplifier and the output signal from the 

10 second amplifier. The output of the subtractor is integrated 
by the integrator, which delivers as output a signal 
representing the flow q(t). Over time, the output signal from 
the integrator, which corresponds to the output signal from 
the calculation resource, converges to a value corresponding 

15 to the instantaneous solution of the differential equation. 

The analogue calculation resource of figure 2 can 
therefore be used to calculate the instantaneous flow and to 
give the direction of flow of a fluid in unsteady motion. 

As illustrated in figure 3, the electronic calculation 

20 resource can also be digital. Figure 3 shows an example of a 
numerical algorithm, programmed in the digital calculation 
resources, designed for the use of a measurement system that 
includes a fluid-velocity measuring device of the narrowing- 
tube type. The algorithm, which relates to the process of 

25 determining instantaneous flow from the pressure difference, 
has 3 stages. The calculation resources include resources that 
are designed to execute these stages. 

A signal is acquired representing the pressure 
difference, every Te seconds. 

30 The first stage 13 consists of determining the initial 

condition allowing the differential equation to be solved. 
Preferably the initial condition chosen will be q 0 = q(t=0) = 
0 when the flow changes direction. 

In the second stage 14, a pressure difference P n is 

35 acquired . 

In a third stage 15, the value q n+i is calculated, 
corresponding to the wanted flow, by discretising the 
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differential equation relating the instantaneous flow to the 
pressure difference, such as by means of the equation : 

q n+ i = q n x (1 - a x | q n | x Te) + P x p n x T e 

5 

At the output of the third stage 15 this gives the value 

q n+ i calculated from q n and Pn . 

By thus reiterating stages 14 and 15, then from the 

differential pressure signal P sampled in accordance with time 
10 base T e , the output q of the digital calculation resource 

supplies a series of discrete values q 0 , qi». q n which converge 

to the instantaneous solution of the differential equation. 

Figure 4 shows an example of results obtained using the 

device of figure 1. The evolution of a pulsed flow in a 
15 conduit is measured by the invention, and compared to the flow 

reconstituted from the speed of the fluid determined by a hot 

wire placed at the centre of the conduit. 

The first graph 17 and the second graph 18 of figure 4 

show changes in the flow in a conduit over time. The dashed 
20 lines show the results obtained with the device of this 

present invention. The solid lines show the results obtained 

with the hot-wire measuring device. 

On the first graph 17, it can be seen that the hot-wire 

measuring device, which is incapable by design of giving the 
25 direction of the flow, always gives a positive flow, even when 

the latter reverses. This is not the case for the system of 

this present invention, which is perfectly suited to reversals 

in the direction of flow. This is due to the fact that the 

measuring system of this present invention takes account of 
30 the direction of the flow. 

On the second graph 18, it can be seen that the device of 

this present invention has a very good resolution in time 

(broad passband) since the pulse rate achieved in the second 

graph of figure 4 is close to 90 Hz. 
35 The pressure measuring system of this present invention 

therefore allows measurement of the flow of a fluid in 

unsteady motion in a conduit. It requires only two static 
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pressure take-offs, and can tolerate any fluctuation of the 
flow, including reversals of direction. This is rendered 
possible by solving the equation relating the flow to the 
pressure difference and used to calculate the flow in real 
5 time, and also gives the direction of the flow. 

The reader will have understood that many variants of 
this present invention are possible. 

For example, the flow-velocity measuring device 1 could 
be a venturi, a diaphragm, or any device that introduces a 

10 pressure drop into the flow. 

In addition, the calculation resource could be of a type 
other than those illustrated in figures 2 and 3, the 
particular feature of this calculation resource being that it 
is designed to solve an equation relating the flow to the 

15 pressure difference, where the said equation is used to 
calculate the flow of the fluid, and takes account of the 
direction of the fluid flow. 

Apart from this, improvements can be made to the system 
and to the process of this present invention. 

20 Flow measurement here has been measured in terms of the 

volume q(t) of an incompressible fluid in unsteady motion in a 
conduit. However, it is possible to calculate the mass flow 
q m , by finding the product of the flow measurement and the 
density p of the fluid. This density is mainly a function of 

25 the temperature of the fluid and the absolute static pressure. 

So by adding a temperature measuring probe to the device 
of figure 1, and using the function p(T) giving the density of 
the fluid according to the temperature in the calculation 
resources, it is possible to calculate the mass flow from the 

30 flow by volume measurement, with the expression: 

q m (t) p(T) * q (t) 

The introduction into the calculation resource of 
35 resources designed to ascertain the density p(T) of the fluid 
from the temperature measurement obtained by means of the 
temperature measuring probe, and resources designed to find 
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the product of the flow by volume q(t) and the density p(T), 
it is possible to calculate the mass flow of an incompressible 
fluid in unsteady motion in a conduit. 

In addition, measuring the absolute static pressure as 
5 well as the temperature enables account to be taken of the 
effects of compressibility of the fluid, and thus to extend 
the applications of this present invention to compressible 
flows (this absolute static pressure measurement can also be 
taken into account when acquiring the density p of an 
10 incompressible fluid so that the mass flow calculation is more 
accurate) . 

Therefore by adding a temperature measuring probe and a 
probe for measuring absolute static pressure to the device of 
figure 1, and by adapting the calculation resource to that it 

15 received as input the pressure difference measured by the 
pressure-difference measuring, resource 2, the temperature 
measurement acquired by the temperature measuring probe and 
the static pressure measurement acquired by the static 
pressure measuring probe, it is possible to calculate the mass 

20 flow of a compressible fluid in unsteady motion in a conduit 
from the equation relating the flow to a pressure difference, 
where the said equation takes account of the direction of the 
fluid flow. 

In addition, the pressure-difference measuring resource 
25 can be other than a differential pressure sensor. For example, 
the two pressure take-offs can be connected to two relative 
pressure measuring probes, so that the relative static 
pressures PI and P2 can then be obtained at the output of the 
probes. The difference between these two relative static 
30 pressures PI and P2 can then be obtained either by means of an 
electronic unit designed to find the difference between the 
two relative static pressures measured, and located between 
the two measuring probes and the calculation resource, or by 
means of the calculation resource 3 if the latter has the 
35 resources to ascertain this difference. Each of the first and 
second pressure take-offs can also be composed of a 
multiplicity of elementary pressure take-offs that are 
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connected together and distributed around the periphery of the 
conduit . 

In this case, the first pressure take-off includes 
several elementary pressure take-offs located in a first plane 
5 perpendicular to the axis of the conduit- The pressure 
measured by means of the first pressure take-off is then equal 
to the mean of the pressures measured by the elementary 
pressure take-offs of the first pressure take-off. Likewise, 
the second pressure take-off includes several elementary 

10 pressure take-offs located in a second plane perpendicular to 
the axis of the conduit, where the first and second planes are 
separate. The pressure measured by means of the second 
pressure take-off is then equal to the mean of the pressures 
measured by the elementary pressure take-offs of the second 

15 pressure take-off. 
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